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Abstract 

Data on Bose-Einstein correlations yield information about the in- 
teraction regions in multiple particle production processes. The con- 
clusions are model dependent. Several popular models are briefly pre- 
sented, compared and discussed. 

PACS 25.75. Gz, 13.65.+i 
Bose-Einstein correlations. 



1 Introduction 

An important reason for studying Bose-Einstein correlations (BEC) in mul- 
tiple particle production processes is that they yield information about the 
interaction regions i.e. about the regions, where the hadrons are produced. 
There is much to learn about such regions. Let us consider two examples. 
In high energy e + e~ annihilations the first stages of the process can be de- 
scribed perturbatively. The two primary leptons merge into a heavy photon, 
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or a Z , this particle decays into a gg pair, the quark and the antiquark 
radiate gluons, the gluons split into qq pairs or into gluons, and so on. Then, 
however, something much more complicated happens: the swarm of parton 
gets converted into hadrons. This is a nonperturbative process, known as 
hadronization, which is poorly understood. One would like to know the size 
and shape of the region where the hadrons are produced, the time interval 
between the first annihilation and the end of the hadronization process, as 
well as the time interval between the beginning and the end of hadronization. 
As a more complicated example let us consider central heavy ion collisions 
at high energy. In this case the initial color fields are so strong that per- 
turbative methods may be unreliable. It is usually assumed that the initial 
partons, whatever their production mechanism, rapidly thermalize. From 
that time on they can be considered as a fluid (a kind of liquid or a kind 
of gas), which expands according to the laws of hydrodynamics or according 
to some simplified version of the Boltzmann equation. Finally the hadrons 
are produced in a hadronization process which may be different from that 
in e + e~ annihilations. Here, besides the hadronization, one is interested in 
the flows of the fluid, in its equation of state and its phase transitions, if 
any. The study of Bose-Einstein correlations among the final hadrons sup- 
plies tentative answers to all such questions. The problem is, however, how 
reliable these answers are? Till a few years ago the mood was optimistic (cf. 
e.g. |X]). The recent results from RHIC, and also some results from LEP, 
however, have been so unexpected and so puzzling (references can be traced 
e.g. from the recent review [2j) that perhaps some important ideas are still 
missing. In the present paper we review some of the most popular ways of 
getting from the experimental data to the physical conclusions. It will be 
seen that much remains to be clarified there. 

The key object in the transition from the experimental data to the phys- 
ical conclusions is usually a single particle density matrix p(pi,p 2 ). There 
are approaches where this is not the case, e.g. the string model (cf. |3] and 
references quoted there), but most models can be formulated in terms of 
this matrix In this section we will present its relation to the data. The 
discussion is grossly simplified (cf. e.g. pQ), but it contains the main ideas. 
The relation between the density matrix and the features of the interaction 
region is the main subject of this paper and is discussed in the following 
sections. Let us stress that the density matrix p(p 1; p 2 ) is an auxiliary con- 
struct, which does not have to coincide with the actual single particle density 
matrix for the mesons of a given type in the final state. 

The diagonal elements of the density matrix are related, as usual, to the 
single particle momentum distribution 
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dN 

^~P(P,P). (1) 

The tilde in most models means equality up to a constant factor. Only 
in the GGLP model it is a more complicated operation. d 3 p may denote 
the infinitesimal volume in momentum space, or the covariant infinitesimal 
volume i.e. the volume in momentum space divided by the corresponding 
energy E(p). The two-particle distribution for identical bosons is 

p(Pi.Pi)p(P2.P2) + Ip(Pi.P2)| 2 - ( 2 ) 



d 3 pid 3 p 2 

For uncorrelated, distinguishable particles the second term on the right-hand 
side would be absent. It results from the symmetrization of the product of 
the single particle density matrices of the two particles and thus reflects the 
BEC. Usually one considers the correlation function obtained by dividing the 
two-particle distribution by the product of the corresponding single particle 
distributions 

rf n r, \ II IP(Pl'P2)! 2 /ox 

P(Pl,Pl)p(P 2) P2) 

In practice the procedure is usually much more complicated than presented 
here, but the strategy is as described: the data are used to get the correlation 
function which is simply related to the single particle density matrix. 



2 GGLP approach 

The first model of BEC in multiple particle production was proposed by the 
Goldhabers Lee and Pais [3] (GGLP). As a realistic description of BEC in 
multiple particle production processes this model is outdated, but it is still 
the best introduction to the subject. GGLP introduced a density p(x) of 
particle sources. They considered identical pions and for definiteness also we 
will call the identical bosons pions, though the analysis can be applied to 
any identical bosons (or even to identical fermions, when the sign between 
the two terms in formula (JBJ) is changed from plus to minus). GGLP made 
the following two assumptions: firstly, that all the pions are produced simul- 
taneously and instantaneously; secondly, that the production is completely 
incoherent - the pions produced in two different space point do not interfere. 
Under these assumptions the single particle density operator is 

p = J d 3 x|x)p(x)(x|. (4) 



3 



The corresponding density matrix reads 



P(Pi,P 2 ) = (Pi|p|Pa) 



d 3 xe~ i( ^p(x.) 



(5) 



where q = p 1 — p 2 . If this were the density matrix which can be obtained 
from formula (j3J), it would be a very nice result. The density of sources 
p(x) could then be obtained unambiguously just by inverting the Fourier 
transform. E.g. for 



Unfortunately, the density matrix (jSJ is untenable as a density matrix pro- 
portional to the true single particle density matrix, because it gives a single 
particle momentum distribution constant in all momentum space. GGLP in- 
troduced, therefore, a projection on the states allowed by energy-momentum 
conservation, i.e. they used the density matrix as calculated here to find 
the momentum distribution for all the particles present in the final state, 
multiplied it by the delta function of energy-momentum conservation and 
integrated over the (covariant) momentum space of all the particles except 
two identical pions for the two particle distribution and over all the momenta 
except one for the single particle distribution. They got results in qualitative 
agreement with experiment. The method, however, was cumbersome - no 
one did calculations for more than six particles in the final state - and was 
in violent disagreement with experiment at the quantitative level 0. 

3 Kopylov and Podgoretskii model 

Very significant progress was obtained by Kopylov and Podgoretskii (cf. [7j 
and references contained there). In their model the density operator in the 
Schrodinger picture is 




(6) 



where R 2 > is a constant, one would obtain 




(7) 
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where 




(9) 



or equivalently 
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where 



xe 



l P x ^(x). 



(11) 



In these formulae x s is a point in space-time, as well as a label which de- 
fines unambiguously a source. It should be related to the position of the 
corresponding source in space-time, but the actual relation is a matter of 
choice. E.g. it could be the point where the source got created, or its mean 
position in space-time. Each source is labelled by a space-time point, but the 
particle produced by a source is not localized in a point. Its wave function 
in space is — x s ) with an additional time dependent phase which is zero 
at t — t s . All these functions are related by shifts in space-time. The wave 
function in momentum space is proportional to A(p), thus all the sources 
yields particles with the same momentum distribution. H is the free parti- 
cle Hamiltonian and we are interested in times t larger than the latest time 
t s , thus the particle evolves freely. 
The corresponding density matrix, 



which can be made to agree with any experimental distribution by a suit- 
able choice of A(p). Thus, there is no obvious need to introduce the energy- 
momentum conservation constraint. Kopylov and Podgoretskii assumed that 
for final states with not too few particles in the interesting momentum region 
the energy-momentum conservation does not affect significantly the one- and 
two-body distributions, and that consequently their matrix ()12j) can be as- 
sumed equal, up to a normalizing factor, to the actual single particle density 
matrix. This assumption makes it easy to get results for high multiplicity ex- 
clusive channels and for high energy inclusive processes. Both were beyond 
the reach of the GGLP method. Besides the satisfactory formula for the 
single particle distribution given above, the model gives for the correlation 
function 





yields the single particle distribution, 




(13) 




(14) 
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Here the factors A(p) cancel and the result is as for point sources. Note, 
however, that because of the four-fold integration this relation cannot be 
inverted to give p(x s ), when the correlation function is known. Moreover, 
a priori the correlation function could depend on the vector q and on the 
vector K = (p 1 + p 2 )/2, while the right hand side depends only on the four- 
vector q. The reason is that the factor A(p) is the same for every source 
and cancels. Physically this means that there is no correlation between the 
momentum of the pion and the position of the source. Since such correlation 
follow from almost every model and since experimentally it is not true, that 
the correlation function depends only on q (cf. e.g. pQ and references quoted 
there) this is a serious weakness of the model. 

4 Yano and Koonin method 

Another attempt to go beyond the GGLP approximations is due to Yano 
and Koonin jS]. In our notation their key formula is 

|p( Pl ,p 2 )| 2 = Re J d^D^pJe^ J d A x 2 D(x 2 ,p 2 )e-^. (15) 
Putting 

D(x,p) = 6(t)p(x) (16) 
one recovers the GGLP model. Putting 

D(x,p) = \A(p)\ 2 p(x), (17) 

which corresponds to the assumptions made in J5|, one reproduces the results 
of Kopylov and Podgoretskii. It is not clear, whether this approach can be 
given sense in the framework of quantum mechanics, when there are position- 
momentum correlations, i.e. when the D function does not factorize into a 
momentum dependent and a space dependent factor. 

5 Covariant current formalism 

The covariant current formalism (cf. 0| and references given there) can be 
derived from the model of Kopylov and Podgoretskii introducing the following 
two generalizations. The label characterizing the source is changed from x s 
to x s ,p s . Thus, the source is characterized by its position in space-time 
and by its four momentum. Correspondingly, the density of sources p(x s ) 
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gets replaced by p(x s ,p s ). Such labels are not subject to the Heisenberg 
uncertainty principle. E.g. a one-dimensional harmonic oscillator may have 
(x) = and (p) = 0. The universal, momentum dependent function A(p) 
gets replaced by the source dependent j(z?r), where m s is the mass andp s the 
momentum of the source. In the Kopylov Podgoretskii model all the sources 
had the same momentum distribution in some overall reference frame, e.g. in 
the center of mass frame of the collision. In the covariant current formalism 
each source has the same momentum distribution, when considered in its rest 
frame. Another way of formulating this model is to replace the incoherent 
sources by incoherent wave packets. Then p s and x s characterize the wave 
packet, and since the waves in the packet describe pions, it is natural to put 
m s = m^. In this model, assuming a fixed mass for all the sources, 

p(pi,P2) = I d 4 x s [ d 3 Ps p(x s ,p s )e^j(^)f(^). (18) 
J J m s m s 

From the point of view of model builders, a nice feature of this approach is 
that one can assume a classical motion of the source, e.g. x s = x s (t s ) and 
Ps = Ps(^s)) and still have a formula which is consistent with quantum me- 
chanics. In the covariant current formalism position-momentum correlations 
are naturally included. 

6 Emission function method 

The method of emission functions is very popular nowadays (cf. e.g. pQ). 
The emission function is built by analogy with the Wigner function ^0], [Hj 
The Wigner function W is related to the single particle density matrix in the 
momentum representation by the formula 

p(Pi, P 2 ) = / d 3 Xe- l( ^W(X, K), (19) 

Where X = (xi +x 2 )/2, and as usual q = p 1 — p 2 and K = (p 1 + p 2 )/2. The 
Wigner function is well-defined as a Fourier transform of the density matrix 
in momentum representation. On the other hand, as a function of X and K 
it may play the role of a phase space distribution. When studying the inter- 
actions regions, it is important to combine information about the space and 
the momentum distributions. Kopylov, Podgoretskii, and followers proposed 
to use the position of the source and the momentum of the particle. These 
two vectors can be measured simultaneously and their joint distribution gives 
an idea about the phase space distribution of particles. The Wigner function 
gives another ersatz phase space distribution. The trouble is, however, that 
the Wigner function refers to a given moment of time, while the hadrons are 



7 



produced during a time interval. The emission function S is supposed to 
improve on that and is related to density matrix by the formula 

P(Pi,P 2 ) = / d A Xe^ x S(X,K). (20) 

The usual strategy is to find from some classical or quasi-classical argument 
a phase space distribution as function of time, and then to interpret it as 
the emission function. Once this is done, one can, by a four-fold integration, 
obtain the density matrix and compare it with experiment. This comparison 
may be used to fix the free parameters of the model. 

Let us note some difficulties of this approach. Given a density matrix 
there is an infinity of very different emission functions, which can reproduce 
it. For instance perfect agreement is obtained for 

S(X,K) = 5(X )W(X,K). (21) 

This, however, corresponds to simultaneous and instant production of all the 
hadrons, which is not a very plausible scenario. When particle production 
is an incoherent sum of production amplitudes at various moments of time, 
the emission function S(X.,t,K) can be related to the Wigner function of 
the particles produced at time t. In the general case, however, the relation of 
the emission function to a Wigner function is hardly visible [TJj. Thus, the 
formula can be used to eliminate wrong models rather than to prove that a 
model is implied by the data. This is not necessarily bad. Eliminating the 
unacceptable values of the parameters of a model, one learns which are the 
good ones. Nevertheless, one must always keep in mind that a completely 
different model can give equally good, or better results. One should also 
keep in mind that X and K are not really position and momentum, but only 
half sums of the corresponding arguments of the density matrix. When the 
contributions of the sources are strongly smeared in coordinate space or in 
momentum space, the difference may be very significant. 
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